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1. Introduction 



Since discovering the first evidence of complexity of phase structure of Higgs- Yukawa 
systems @, @ much has been done to understand them better (for a review and list 
of references see [|J, 0). There were studied phase structure and continuum limits of such 
systems with different lattice formulations and symmetry groups; proposed hypotheses on 
possible non-trivial fixed points and checked some of them. There were used numerical, 
analytical, and combined methods. Much has been understood about this, much is still to 
be answered. In particular, such questions as: which features of the phase diagrams are 
inherent properties of the Higgs- Yukawa model, and which are lattice artifacts; whether 
non-trivial fixed points exist at intermediate values of the Yukawa coupling; if so, what 
does distinguish such points on the phase diagrams and whether this depends on group of 
symmetry of the system, have not yet definite answers. 

This paper, not claiming to solve definitely those unanswered problems, rather pursues 
the following aims: to give a simple analytical technique for examining on the same basis 
the phase structure of a wide class of lattice formulations of the Higgs- Yukawa systems with 
various symmetry groups; to compare the phase diagrams for some of such formulations 
for Zii U(l), and SU(2) symmetrical systems; to give additional arguments pro et contra 
those previous results which seem to be not well-established, and to draw one's attention 
to interesting points which were left being not understood well. 

We shall use variational mean field approximation and analyse the systems with ra- 
dially frozen Higgs fields and chirally invariant lattice fermion actions possessing real 
fermionic determinants. These conditions enables us to calculate contributions of fermionic 
determinant into mean field free energy — the main problem of the method, in a ladder 
approximation ||, in other words, to sum up the contributions of leading and some of 
those of next-to-leading orders in inverse space-time dimension 1/D knowing only form of 
the free fermion propagator in the momentum space. Compared with other mean field cal- 



culations M, M, M, [1C], [11] this gives us two advantages: to analyse the phase diagrams 



of the systems with a wide class of lattice fermion actions, and, in fact, for any values of 
the Yukawa coupling, including the most interesting region of intermediate ones. On the 
other hand the above conditions do not allow us to discuss interesting features revealed in 
a model with non-frozen Higgs field |TI| and to compare our results with those obtained 
with staggered fermions @, |13[ (except formulations of |§, HlPf , [jlTf with a number of 
staggered fermions is a multiple of 2 15 / 2 ), and with results for the Wilson fermions 0, [II 
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By virtue of the Nielsen-Ninomiya theorem [|15], chirally invariant lattice fermion 
actions (providing they are bilinear in fermion fields and lattice translation invariant, the 
conditions we shall respect) must be either non-Hermitean, or non-local, otherwise they 
involve an equal number of the left-handed and right-handed Weyl fermions. Using our 
technique we cannot take non-Hermitean actions, being able to consider naive, non-local, 
and mirror fermion actions. In all the cases the fermions coupled to the Higgs fields in the 
same local way. 

The outline and the results of the paper are as follows. 

The systems under consideration are defined in Sect. 2. We use their usual lattice 
parameterization in terms of scalar hopping parameter k and Yukawa coupling y. In Sect. 
3 we describe the method and approximations, and obtain the explicit formulae for second 
order phase transition lines n cr (y) for the symmetry groups Z 2 , U(l), and SU(2). 

In Sect. 4 we apply these formulae to the systems with naive, SLAC, a maximally 
non-local (the Weyl), and a mirror fermion actions. We find that the phase diagram for 
all the systems show a kind of universality at k > 0, while at k < they crucially depend 
on both fermion action and symmetry group. In particular, ferrimagnetic phase does not 
appear in any systems, while appearing in £7(1) ones with the Weyl and mirror fermion 
actions, and in all SU (2) systems. 

In Sect. 5 we make mean field estimations of fermionic propagators and conden- 
sates along the ferromagnetic-paramagnetic phase transition lines. Both show different 
behaviour in weak and strong coupling regimes in all the systems. This enables us to 
speculate on locations of possible non-trivial fixed points. We conclude that the most 
interesting ones can exist in Z2 system with SLAC action and in U(l) systems with naive 
and SLAC actions. We also note impossibility of defining a chiral theory in paramagnetic 
phase of mirror fermion model at strong coupling. 

Sect. 6 is a summary and discussion. 

2. The system 

The system is defined on a hypercubic D-dimensional (D is even) lattice A with sites 
numbered by n = (ni, ...,ud), —N/2 + 1 < < N/2 (N is even, and eventually tends 
to infinity) and with lattice spacing a = 1; ft is the unit vector along a lattice link in the 
positive /x-direction. Dynamical variables are Dirac fermion fields ift n , ip n , and scalar fields 
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$ n G G, where group G = Z 2 , U(l), or 577(2), so that = 1- We use the following 

representations for the group elements: 

p 

$ = T</> = ^T a (/) a , (2.1) 

a=0 

with real </> a such that Y7a=o^n) 2 = 1- Hence we have 

p = 0, T = l, Z 2 , 

p = l, T = (l,i), (2.2) 

p = 3, T = (1, zr 1 , z'r 2 , ir 3 ), SU(2). 
We imply antiperiodic boundary conditions for the fermion and periodic for the scalar 
fields. 

The system is defined by functional integrals 

Z[J] = / ]_]_d$„# n # n e (2.3) 

with action 

A[$,lP,tf\=-2K K0n+» + J2^rn[0mn + y(PL$l + PR$ m )S m n}lpn, (2.4) 

n,ii,a m,n 

where <i$ n is the Haar measure on G; k G (—00,00) is hopping parameter; y is the 
Yukawa coupling which without loss of generality will be considered non-negative; Pl,r = 
(1 ± 7£>+i)/2 are chiral projecting operators; is a lattice Dirac operator determining the 
form of the fermion action (the systems with mirror fermions which we consider in Sect. 
4.4 is reduced to this form). 

We shall consider actions with operators satisfying the properties 

0mn 0nm i 

f zp(m — n) ^-^ 
(Pmn = / e (2.5) 

L l (p) = l h ip) , L »(-p) = ~ L M> 

where J = J d D p/(2n) D , Pfl G ( — 7T, 7r), and we use the Hermitean 7-matrices: [7^, 7^]+ = 
2^. 

Action (|2.4| ) is invariant under G x G global chiral transformations 

-> (h L P L + h R P R )ip n , 
^ n ^^ n (P R h{+P L h R ), (2.6) 
$ n -> /i L $ n /ijj, 

with Iil r £ G. 



3. Method and approximations 

To analyze the phase diagrams of the system we use the variational mean field ap- 
proximation |16|] which becomes applicable to ( |2.3| ) after integrating out the fermions 

-W[J] 

(3.1) 



Z[J] = e 



= /n^ 



2«E re ,,, fl r n r n+ a + Indet [0 + y*] + E n ,a « 



where $ = (Pl& +Pr^)- Then for free energy of the system F = W[0] the method yields 
inequality 



F < Fmf = inf 



J>(M +J2 h n(<t>n)h 
n n,a 

-<2« 0; +A + lndet[^ + j/f]), 



n,fi,a 



(3.2) 



1/2 

where fo^ is the mean field with radial component h n = [E a (^™) 2 ] > 

Ea 



-In / d$„e 



= - In cosh ^ = - l -h 2 n + ±-hi + 0(h 6 n ] 



ha I (h 7 



= - In— Ji(/i r 



lh 2 n + ±hi + 0(h 6 n ), U(l), 

= -l h n+^ h n + °(h 6 n), SU{2), 



(3.3) 



is the test free energy per lattice site [I v (z) is the modified Bessel functions], and 



Hd$ n O[<j>}e 



En,a KK 



(3.4) 



The inequality (|3.2|) is expected to tend to equality in the limit of D — ► oo fl6[j . So, 
we can get some idea of the system at D = 4 studying Fmf- From (|3.3| ) it immediately 
follows that 







u(h r 



Oh*" ' " 
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K + Oih^), z 2 , 
lh a n + 0(hl), U(l) 



(3.5) 



= -h a n + 0(h 3 n ), SU(2); 



= h m h n + 5 mn (l - hi) +0(h*), Z 2 , 

= ^^n + ^mn[5 a6 (l-^ 2 J-^>y+C»(^), 17(1), 

= ^KJl + \5 mn [5 ab (l - ±h 2 n ) - + O^}, 517(2); 

and so on. Therefore the main problem is calculation of expectation value (lndet [0+yQ])h- 
In our case it has the following representations 

(lndet [0 + y$\) h 

1 00 1 

= lndet$]--£-y 2 * £ tr [(^ 0~] 3 • • -^X^ • • • ^) h ] 

n=l ii,... ,i 2 n ,„ „x 

= c 9 2 D / 2 iV D lny 1 ' ' 

1 °° 1 1 

n=l H,---,«2tj 

where tr stands for the trace over spinorial indices, as well as for group ones for SU(2) in 
which case factor c g = 2, otherwise c g = 1, and it has been taken into account that the 



trace of an odd number of 7-matrices vanishes; correlators of the Higgs fields in ( |3.7|) in 
terms of their real components read as 

= E T^ t T^---T a -(^c 2 2 ---c;)^ ( 3 - 8 ) 

We consider F^f on translation invariant ansatz for 

h a n = h a + e n h a st , e n = (-l)^ n ». (3.9) 
Then the mean field equations are reduced to 

■777 F MF = 0, — — Fmf = 0. (3.10) 

oh oh st 

In view of (|3?5| ), the solutions of Eqs. ( |3.10|) , h* ("magnetization") and h* st ("staggered 
magnetization"), in fact are the order parameters distinguishing the ferromagnetic (FM): 
h* 7^ 0, h* st = 0; antiferromagnetic (AM): h* = 0, h* t 7^ 0; paramagnetic (PM): both are 
zero, and ferrimagnetic (FI): both are nonzero, phases in the system. 
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Further simplification comes from the observation (see, for example H), that as the 
values h = 0, h st = are always solutions of Eqs. ( |3.10| ), and, therefore, the second order 
phase transition lines are determined by equation 



Q2 Q2 Q2 



mf or 9 -Fmf — or or — Fmf = (3.11) 



at h — 0, h* t and h* , h s t = 0, in order to find the critical lines not far from the PM phase, 
it is sufficient to know Fmf to terms of order of h^. 

If the problem could be solved exactly both of two representations ( |3.7| ) of the fermionic 
determinant would yield the same answer. But correlations of 0'S at coinciding arguments 
make the problem unsolvable exactly. Indeed, the contributions of order of to (|3.7| ) 
come not only from terms oc u' 2 [first terms in Eq. ( 3.6|) 1, but also from terms of any 
orders in u" [terms proportional to S mn in (|3.6| )1, as well as from higher correlators. Some 
of such contributions shown schematically in Fig. 1. Therefore, we are forced to use 



some approximations, and, in particular, to use two representations of (|3.7| ) separately for 
"weak" and "strong" coupling regimes of y, though the exact meaning of this can only be 
clear a posteriori. 

Our approximation involves summing up all diagrams of Fig.l (a) (proper lad- 
der diagrams) and (6) (crossed ladder diagrams), so we call it ladder approximation. 
Diagrams (a) and (b) correspond to contributions to (lndet[$ + y$])h °f the form 

u 'h u i n +i ( u ") 2n ~ 25 i2i2 n ---Ki n+ 2i and i u ") 2n S hin+i 5 hi 2n ■ • ■ S i n in+^ respectively. 

Terms oc hh st do not contribute to Fmf (due to momenta conservation). Hence, Eq. 
( [3. 11|) in neighbourhood of the PM phase turns to pair independent equations 

d 2 d 2 

F MF\ h= u =r ,= 0, -^rTF MF \ = 0. (3.12) 



dh 2 < h=h s t=o 9/,r*'ik=k„=o 



Then, using properties ( ^.5| ) of the Dirac operator, formulae ( |3.6| ), and the properties of 
operators T: 

v 

T at T b T aj = 1, Z 2 , 



a=0 

= 0, U(l) 



(3.13) 



-2T b \ tr T at T 6 = 25 ab , SU(2), 



we find that the second order phase transition lines for the system ( |2.3| ), ( |2.4| ) in our 
approximation are determined by the expressions: 
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fMW (y) = _L| 1 _ 2 ^ 



y 2 G w (0) 
l + y 2 G w (0) 



ff y 2 G w (g) 4rW 2 



i M{w) (y) = -^{i-i D/2 



FM(S ){y)= _l_^_ 2D/ 2 



y 2 G w (n) 
l + y 2 G w (n) 

y 2 G w (g) 

(l + y 2 G^(q)) 2 
G s (0) 



+ y 4 G w "( q ) 



y 2 + G s (0) 

y 2 G s (q) 



K 



AM(S) 



(y) 



__L/i_ 2 °/ 2 



4D { 



(y 2 + GS(q)) 2 y 
G s (n) 



+ ~,G S \q) 



y 2 + G s {ir) 



2 G s (q) 



(y 2 + GS(q)) 2 y 



U(l): 



4r M{w \y) = ^V-z D,2 trG w {v) 



*t M(s \y) = -^ 



1 

i 

2? 



l_2°/ 2 -^G 5 (0)], 



1 _ 2^/ 2 JL G s {n) 



SU(2) : 

: FM{W)^ 



K, 



1 

D 



1 - 2 D ' 2 



K 



AM(W) 



(») = — 1"2 D/2 



FM(S) 



(y) 



D 



y 2 G w (Q) 

-2/^(0) 

2y 2 G ,w/ (g) 
(2-y2 G w (g)) 2 

y 2 G w (7r) 
-y 2 G w (7r) 
2y 2 G w (g) 

G s (0) 



1 

I 1 



7^ 2 ( ) 



(3.16) 



AM(S) 



(y) 



i 



1 - 2 D ' 2 



2y 2 -G s (0) 

2y 2 G s (q) 
(2y 2 -GS(q)) 2 

G s (ti) 



(?) 



where 



G W (q) 
G S (Q) - 



2y 2 -G s (n) 

2y 2 G s (q) _ 
(2y2_ G S( g ))2 

L(p)L(p + q) 
L 2 (p)L 2 (p + q) 

L(p)L(p + q), 



(3.17) 



(27T) 



^ G (-7T,7r]. 



In (|3.14|) - ( |3.16|) functions n FM (y) describe the critical lines between the FM and 
PM phases, and k am (y) those between the AM and PM phases, while superscripts W 
and S mean strong and weak coupling regimes, respectively. The contributions to ( |3.14| ) 
- ( p.!6|) which are proportional to G(0) or G(ir) come from the diagrams of Fig. 1(a), 
while the integral terms from those of Fig. 1(6). The second terms in the integrands 
take into account that the first crossed ladder diagram enters with a factor 1, rather than 
2 (in preprint version of || and in |TJJ such terms in the equations for Z<i system were 
erroneously missed). In the case of U(l) due to the specific property of operator T ( |3.13| ) 
both ladder and crossed ladder diagrams vanish, and all the contribution of the fermionic 
determinant is due to the first diagram of Fig. 1(a). In the case of nf fermions couple to 
the Higgs fields with the same y, factor nf2 D / 2 appears instead of 2 D / 2 in Eqs. ( gTj ) - 
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The diagrams of Fig. 1(a) are generalization of "double chain" diagrams of Refs. ||, 
||1Q| to any configurations of the same topology. These ladder diagrams coincide with the 
double chain ones only in the case of strict locality of the Dirac operator (or its inverse 
operator). In fact, this is the case only for the strong coupling regime for the systems with 
naive action (see Sect. 4.1). Crossed ladder diagrams of Fig. 1(6) are generalization of 
the double chain diagrams with coinciding ends, which have not been taken into account 
in previous calculations. Although the crossed ladder diagrams are of 0(D~ 1 ) compared 
with the proper ones, and can be neglected at "weak" as well as at "strong" couplings, 
they can become dominating at "intermediate" ones, when y is close to the singular points 
of the integrands in (|3.14j ) and ( |3.16|) . 

are basically determined by four constants G w (0), G w (ir), G s (0), and G s (n) 1 which, 
in their turn, are determined by the form of the free fermion propagators. It is these 
constants that determine the singular points of the expressions for Zi and SU(2), thereby 
determining domains of the weak and strong coupling regimes in these cases. As G(0) > 0, 
G(tt) < 0, these domains are: 

Z 2 : 

y<y w = \G w (7r)\-K y>y 3 = \G s (7r)\^, (3.18) 

SU(2) : 

y<y w = [^(0)p, y>y s =[±G s (0)]K (3.19) 

There can be the following possibilities: 

(i) y w < y s . Although we have no analytical expressions describing the system in 
the region y w < y < y s ', we can get, in fact, complete picture of the phase diagrams if the 
integral terms diverge at the points y w and y s and no FI phases appear. These are the 
cases of Z 2 systems with naive and mirror fermions [Figs 2(a), (d)]. 

(ii) y w > y s . The domains of the weak and strong coupling regimes are overlapped, 
so we can continue the lines k w (y) and K S (y) until they intersect each other. In these 
cases we know the phase diagrams in fact at any y [an example is U(l) system with naive 
action, Fig. 3(a)] (see, however, remark at the end of Sect. 6). 

(Hi) Lines k fm intersect lines k am forming the FI phases before cases (i) or (ii) 
are realized. In such a case formulae ( |3.14j ) - (|3.16| ) describe the phase diagrams only 
in a neighbourhood of the PM phases (examples are all SU(2) systems, Figs. 4). The 
important fact is that the lines together with their first derivatives are smooth in the 
points of intersection. This follows from Eqs. ( |3.12| ). 
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For U(l) systems we continue lines k w (y) and k (y) until either case {%%) or (Hi) is 
realized. In this case, however, such a procedure need to be justified. The natural and 
simple way is checking that the intersections of the curves occur at the points at which 
quantities (y 2 /2)G w and G s /(2y 2 ) are obviously less than 1. 

The proper ladder diagrams give the contributions beginning from the order 0(?/ ±2 ), 
the crossed ladder ones from 0(y ±4 ) (note that f G(q) = 0). Contributions of other 
diagrams [Fig. 1(c)] come into play in higher orders in ?/ ±2 , at least from the order of y ±6 . 
So, we expect that these contributions are non-singular and numerically suppressed. 



4. Phase diagrams 



In this section we apply formulae ( |3. 14j) - ( 3.16Q to four-dimensional systems with 



different chirally invariant formulations of fermions on a lattice. All the phase diagrams 
are shown for n/ = 1, unless other is indicated. 

4-1. Naive fermion action 

In this case operator is local 

, 2 (4.1) 
L M (p) = sinp M , 

and therefore produces species doubling. The system is invariant under transformations: 
(ip,ip) n -> exp(ze n 7r/4)('0, ip) k —k, y — —iy M, that, in particular, 

results in G(n) = —G(0). Numerically one has 

G w (0) = -G w (n) w 0.620, G s (0) = -G s (n) = 2. (4.2) 

Phase diagrams for Z 2 , U(l), and SU(2) at rif = 2 are shown in Figs. 2(a), 3(a), and 
4(a), respectively. Some comments follow. 

Zi'- Fig. 2(a). Domains of the weak and strong coupling regimes are determined by 

y w w 1.27, y s w 1.41. (4.3) 

The lines nf r M (y) and k^ 4 (y) do not intersect each other, so no FI phase appears, the case 
(i) being realized. The lines form two disconnected domains with PM phases, as well as 
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with AM phases; the funnel with FM phase is formed by logarithmic dropping of the critical 
lines near the points y w and y s : K^ r M{W) oc ln[l + y 2 G w {>*)}, K f r M(S) oc ln[l + G s (n)/y 2 ], 
and therefore extends up to k — > — oo. 

U(l), rif = 2: Fig. 3(a). The lines form three connected domains with FM, PM, 
and AM phases; no FI phase appears. These features are hold for any n/. Both curves 
intersect each others (point A and its counterpart on PM-AM line) at y m 1.34, where 
the construction is justified. The result agrees with previous analytical ones |J, |JTT|| , 



but disagrees with numerical results of Ref. , where an evidence of FI phase has been 
revealed. This fact, as well as point A which is a candidate for non-trivial fixed point, is 
discussed in Sect. 5 and 6. 

SU(2) 1 nf = 2: Fig. 4(a). Domains of the weak and strong coupling regimes are 
overlapped: 

V W « 1-80, y s = l, (4.4) 

the case (Hi) being realized. The FI phase appears from rif = 1 extending with rif. In the 
domain of applicability of our formulae (solid lines) the phase diagram is in a quantitative 
agreement with the Monte Carlo results of Ref. . 

4-2. SLAC fermion action 

This action is defined by operator of the form |18 



X^( -1 ^ +1 7 (dm+lp. n ~ Sm-lp, n), 



i>o 1 (4.5) 

L^(p)=P^, PfM e (-7T,7r). 
It represents an action with a moderate non-locality as mn drops with distance like \m — 
n\ . In this case we have 

G w (0) w 0.109, G w (n) w -0.0544, G s (0) = ^n 2 , G s (n) = -\ 2 . (4.6) 

Phase diagrams are shown in Figs. 2(6), 3(6), and 4(6). 

Z2: Fig. 2(6). Domains of the weak and strong coupling regimes are overlapped: 

^^4.29, y s ^2.57. (4.7) 

There are three connected domains with FM, PM, and AM phases, and no FI phase. These 
features are hold for any rif. 

U(l): Fig. 3(6). The phase diagram for n/ = 1 looks like in the case of Z2. But for 
rif > 2 the FI phase appears. Points of intersection of the critical lines located at y ~ 3.31, 
where the construction is justified too. 

SU(2): Fig. 4(6). Domains of the weak and strong coupling regimes coincide with 
(for y w within calculational errors). FI phase appears for any rif > 1. 
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4-3. Weyl fermion action 

This action is denned by the finite dimensional approximation of functional integrals 
for Weyl quantization |T3]. In this case we have 



5>MH(-l) l+1 2(*m+lA»-* 



m — lfi n) j 



l>0 



, (4-8) 
L M (p) = 2tan-p^. 

This is a maximally non-local action in the sense that mn does not drop at all with 
increasing \ m—n\. But a remarkable fact is that this action can be transform to a local form 
if we introduce variables ip defined on the centres of D-cells of the lattice, i.e. on the sites 
of the dual lattice, leaving fields i/j being defined on sites of the original one. Then change 
of variables, which in momentum space looks like ifi p = F(p)ipp, with F(p) = [7 cos \pp, 
leads to a local action with L^(p) = 2 sin \p^ Yl^^ cos \Vv- Although now L^(p) has 
additional zeroes at the Brillouin zone boundary, the system of course is not changed: 
contributions of the additional species to the partition function are canceled by Jacobian 
coming from the change of variables, while their coupling to the Higgs field is suppressed 
by the factor F(p) (in this point this looks similar to the Zaragoza proposal ||2U[|). 

The more non-locality of the action, the less and the grater \G S \. In this 

extremal case we have (at N — > oo) 

G w (0) « 0.0450, G w (tt) « -0.00739, G s (0) -> oo, G s (n) = -16. (4.9) 

The divergence of G s (0) means that in this case terms of the strong coupling expansion 
diverge, and cannot be summed up into the finite expression in the ladder approximation 
(in [17 1, because of missing the second term in the integrand of ( |3.14| ), the wrong conclusion 
has been made on this point). Thus, we can analyse the phase diagram only in the weak 
coupling regime; the results are in Figs. 2(c), 3(c), and 4(c). 

Z2: Fig. 2(c). Formally domains of the weak and strong coupling regimes are over- 
lapped: 

y w &11.6, y S = 4. (4.10) 

Despite a tendency, no FI phase appears. 

£7(1): Fig. 3(c). The FI phase appears at y ~ 5.16, where our condition is satisfied. 
577(2): Fig. 4(c). In this case formally we have 

y w w 6.67, y s -»■ 00. (4.11) 

The phase diagram is similar to that of £7(1) case, however, in view of ( |4.11| ) it looks 
plausible that at y > y w only the FM phase exists. 
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4-4- Mirror fermion action 



We consider the simplest variant of the mirror fermion action [2~I| with zero (bare) 
mixing parameter between fermion field ip and its mirror counterpart x an d with only X 
coupled to the Higgs field 

A = E 9m (Cn ^n + W mnX n)+Xm ^n + W mn ^ n )] 

(4.12) 

where $ N is Dirac operator for naive fermions ( f4.1| ), while is the Wilson operator 

W mn = — — (<5 m +/i n + 0~m-jl n — 2 8 mn ). (4-13) 

The action has the mirror symmetry 

ljj n -> (/l L P L + h R P R )^ n , ^ n -> ^ n (P«4 + 

Xn -> (/iflPi + h L P R ) Xn , Xn ~> X«(-Pr>4i + ^4)> ( 4 ' 14 ) 
$ n -> h L $ n h R . 

We can apply our formulae to this system after integrating out ift. Then, as ip does not 
couple to <p and therefore its determinant is an irrelevant constant, we come to effective 
non-local action in terms of fields x an d 4> °f the form of (|2.3| ) - ( |2.5| ) with 



iN W(0 N )~ 1 W] 



mm 



Lp(p) = sinp^ 



1 (E^ 1 -cosp^))' 

• 2 



(4.15) 



This non- locality is of a new type compared with two preceding cases as L n (p) now 
involves also all p v with v ^ \x. Now we have 

G w (0) w 0.0259, G w (7r) w -0.00734, G s (0) w 348, G s (7r) w -159. (4.16) 

Phase diagrams are shown in Figs. 2(d), 3(d), and 4(d). 

Z2: Fig. 2(d). Domains of the weak and strong coupling regimes are not overlapped: 

y w ss 11.7, y s ss 12.6, (4.17) 
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so this case incorporates features of both local and non-local actions. The phase diagram, 
except the scale of y, has the same features as that for naive fermion action. 

U(l), rif = 2: Fig. 3(d). Domain with FI phase appears from rif = 1, expanding with 
rif. For rif = 2 it appears at y ~ 5.19 and y ~ 9.72, so that according to our criterion we 
can trust the picture. The phase diagram qualitatively agrees with the result of Ref. [p2 | 



(there the phase diagrams examined in a region of parameter space that is different from 
ours) . 

SU(2): Fig. 4(d). The weak and strong coupling domains are determined by 

y w « 8.79, y s « 13.2. (4.18) 

As it happened in all other SU (2) systems, FI phase appears. 

To show the relative role of the proper and crossed ladder diagrams, we display in 
Figs. 2 (by thin dashed lines) contributions of the proper ones. This demonstrates that 
crossed ladder diagrams play the important role only when the case (i) is realized, being 
in fact negligible in other cases. We expect that contributions of diagrams that we did not 
take into account are actually invisible, at least in the cases (ii) and (iii). 



5. Fermion correlators 

In order to learn more about structure of various phases of Figs. 2 - 4 we now make 
a mean field estimation of fermion correlators (ip m ift n ) for the above systems. We mainly 
concentrate on an neighbourhood of the PM-FM critical lines, as it is this domain that is 
expected to be the most interesting for the continuum physics. 

We shall evaluate the quantity 

(tilMF^^ + t/^^,, (5.1) 

where expectation value in the r.h.s. is defined in ( |3.4|) , while h* is a solution of the mean 
field equation ( |3.1U| ). Such an expression appears very naturally as a variational mean field 
approximation for the correlator (see ||), though there is no strict relation similar to Eq. 
( p.2| ) for free energy. By definition this is a quenched estimation. 
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We have the following representations for (|5.1| ): 

oo 



-J ("'•- 1 

.. oo 



oo 1 

'$t \ 5 + Vf-iV— — V $ ■ $ ! ■■■$■ 



We shall use these two representations for weak and strong coupling regimes which have 
been determined by our preceding considerations. 

Consider first (V'mV'n) mf in FM phase in the approximation of uncorrelated Higgs 
fields. Choose their expectation values to be real: (4>n)h* = 8 a0 ((f)), so that 

($! $ 2 ■ • • $,> h . = ($i $j ■ • ■ $+> h . = (<f>y. (5.3) 



Then, from (|5.2|) it follows: 



T (5-4) 
This reproduces the well-known result, which has been obtained by various methods (for 



the first references see [p3 [, , ||13|1 )) that the behavior of the fermion masses with (</>) (and, 
therefore, that of renormalized y) is completely different in the weak and strong coupling 
regimes. While in the weak regime one has the usual perturbative Higgs mechanism (with 
the Gaussian fixed point for y), the fermion masses at strong coupling do not vanish on 
the critical lines and in PM phase. In this approximation they tend to infinity, that is, the 
fermions decouple there. We shall refer hereafter to PM (FM) phases at weak and strong 
coupling regimes as PM(W) and PM(S) [FM(W) and FM(S)], respectively. 

In [17j it was contemplated a possibility to use this feature of the Higgs- Yukawa 



systems for defining a chiral theory in PM(S) phase of the mirror fermion model ( ^.12|) . 
The idea was that the mirror fermions x decouple there leaving behind massless chirally 
invariant fermions i/j. Using the above approximation, however, it is easy to show that 
though all goes in such a way, the goal cannot be reached: ip turn to naive fermions. 

Although approximation ( |5.3| ) yields correct qualitative picture, it is too rough for a 
quantitative analysis: it has been shown by numerical calculations ffl, [EI, M, that the 
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fermion masses at strong coupling indeed increase with decreasing (</>), but remain finite 
on the critical line, rather than tend to infinity. 

To proceed further we make the mean field estimation of the fermion condensate 

l^^)MF = N- D J2@nlPn)MF (5.5) 

n 

along the FM-PM critical lines where it has the form: 

®^)mf = ~2 D/2 C(y) (<P) + 0((</>> 2 ). (5.6) 

This allows us to use the ladder approximation, that is, to sum up the contributions to 
) of the diagrams of Fig. 5 (a) . Then, from (|3.5| ) , ( |3^6| ) and ( |5.2|) we find 

yG w (o) 



C&(y) = 



l + y 2 G w (0) 



y 



(5.7) 



y 2 + G s (0) 



17(1): 



SU{2) 



C^(y) = yG w (0), 

y 



(5.8) 



2-y 2 G w (0) 

2y 



(5.9) 



2y 2 -G s (0)' 

In the case of U(l) only the first diagram of Fig. 5(a) gives a contribution to ( |5.8|) : the 
situation is very similar to that for critical lines (|3.15| ) (see also 

condensates are smooth functions along both PM(W)-FM(W) and PM(S)-FM(S) crit- 
ical lines. In the cases when FI phase appears it prevents us to follow far beyond the 
points of intersection of PM-FM and PM-AM lines. The important fact, however, is that 
in a neighbourhood of those points condensates remain smooth functions of both their 
arguments y and (<p). This is due to the fact that the staggered magnetization gives no 
contributions of order O (((/))) to (i/ji/j)mf (because of momentum conservation). 

In the cases when the PM-FM line is continuous and no FI phase appear, the function 
C(y) can be discontinuous at the points of intersection of FM(W) and FM(S) lines. We 
have only three such systems: Z 2 with SLAC fermions, and U(l) with naive and SLAC 
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(rif = 1) fermions [points A in Figs. 2(6) and 3(a), (6), respectively]. Figs. 6 and 7(a), (6) 
show that this is indeed the case. This means that the condensate, being zero in both 
PM phases and on the whole critical line, is discontinuous in FM phase. The discontinuity 
of the fermion condensate is an evidence of the first order phase transition. Indeed, the 
condensate can be defined as the first order derivative of the free energy in respect to an 
infinitesimal fermion mass (see also 0). As the condensate is an order parameter of the 
systems, points A in Figs. 2(6) and 3(a), (b) look like tricritical points in which the first 
order phase transition turns to the second order one. This allows us to identify the points 
A with point A in Fig. 2 of Ref. |l] , thereby considering them as candidates for non-trivial 
fixed points. 



6. Summary and discussion 

We derived the explicit formulae [Eqs. (|3.14|) - ( |3.16| )1 describing phase diagrams of a 
wide class of Z%, U (1), and 577 (2) Higgs- Yukawa systems, and applied them to the systems 
with naive, SLAC, the Weyl, and mirror fermion actions (Figs. 2-4). The phase diagrams 
turned out to be very different for different symmetry groups and fermion actions at k < 0, 
being of the same form at k > 0. 

The difference between the phase diagrams for different formulation of lattice fermions 
shows a lack of universality in the systems at k < 0, and can be interpreted as a lattice 
artifact. It is well known that in this region the sufficient condition of reflection positivity 
is not satisfied, and, therefore, relevance of the systems to well-defined quantum field 
theories is under the question. If, however, in that region physical positivity is fulfilled, 
the systems are still interesting from the point of view of continuum physics. Analysis 
of scalar propagators at k < in 577(2) system with naive fermions 0, |23|] gives some 
evidence that this is indeed the case. 

All the points which are candidates for non-trivial fixed points lies in the region of 
negative k. There are two types of such points: those where PM critical lines intersect AM 
lines forming FI phase, and the points where we can expect a phase transition separating 
weak and strong coupling regimes. The points of the first type present in all SU (2) systems, 
Figs. 4, and in some of U(l) systems, Figs. 3(c), (d). This case, however, seems to be 
excluded: the numerical investigations of such points in 577(2) system with naive fermions 
0,0, gave no evidence of non-trivial behaviour of the system. Our analysis speaks 
in favour that too: we applied fermionic condensate calculated along the PM-FM critical 
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lines [Eqs. (|5.5|) - (|5.9|) ] as an order parameter sensitive to first order phase transitions, 
and showed that the condensate, as well as both critical lines, does not feel these points, 
remaining smooth functions of their arguments. 

Among the points of the second type only those are relevant to continuum physics 
that border on PM phase: at other points either (0) or ((j) s t) is not zero. Therefore, the 
most interesting are systems which have continuous PM-FM critical lines and no FI phase. 
In our examples these are Z2 system with SLAC action and U(l) systems with naive and 
SLAC (rif = 1) actions. Indeed, in FM phase near such points the fermion condensate is 
discontinuous, Figs. 6 and 7(a), (6), that is an evidence of the first order phase transition 
separating the FM(W) and FM(S) phases. We therefore identify the points A in Figs. 1(6), 
and 2 (a), (5) with point A in Fig. 2 of Ref. ]lfl. If these points are really non-Gaussian, than 
non-trivial continuum theories with two relevant parameters, Higgs and fermion masses, 
can be defined approaching the points from FM(W) phase. 

The SU(2) systems we considered have no such points. The latter, however, can exist 
in SU(2) systems with other formulations of lattice fermions, in particular, in those with 
staggered fermions coupled to the Higgs fields in a local way, when number of the fermions 
is not too big [|IU|, [TJJ- 

We expect that disagreement of the phase diagram for U(l) system with naive fermion 
action, Fig. 3(a), with results of Ref. @, where an evidence of FI phase has been found, is 
due to finite lattice effects in the numerical calculations. Our results for SU(2) system with 
naive fermions are in quantitative agreement with numerical calculations of Ref. , and 
results for £7(1) system with mirror fermions are in qualitative agreement with those of Ref. 
|22Tl (we cannot compare them directly, because the latter have been obtained in different 
region of parameter spase of the model). In both cases the appearance of FI phases has 
very clear reasons. Therefore we do not see any reasons why our approximation could 
fail in this case. The more so, that in this case it is most stable against 1/D corrections 
(cf. 0). Otherwise, there must be some underlying physics which by unknown reasons is 
not taken into account by our approximation, and which worth further investigation. The 
only change of the above picture which we cannot exclude is that the FM-PM (perhaps 
together with AM-PM) critical line is actually discontinuous, being teared up by the first 
order phase transition line at some value of y not far from y A « 1.34. Then the point A 
splits into two first order phase transition points where critical lines ttcr M ^ W ^ and Kcr M ^ S ^ 



end up. In this case, in view of the results of Ref. [12], one can hardly expect the existence 
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of non-trivial fixed points in the system. Examination of this issue, however, requires 
another technique. 

Therefore, it would be very interesting to repeat thorough numerical investigation of 
the U(l) system with naive fermion action near the point (y A , k a ) ~ (1-34, —0.43) (for 
nf = 2). Another interesting issue is to trace the evolution of the phase diagrams of such 
system at finite scalar self-coupling A |T5[ with increasing A. 
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Figure captions 



Fig. 1. Diagrams contributed to the expectation value of the fermion determinant 
(|3.7|) to the order h^', (a) is the ladder, (6) the crossed ladder diagrams. Solid lines stand 
for or 0~ x , solid circles for u', dashed lines for u" . 

Fig. 2. Phase diagrams of Z<i systems with (a) naive, (b) SLAC, (c) the Weyl, (d) 
mirror fermion actions. Dashed lines show the contribution of only diagram of Fig. 1(a). 
Point A in (6) is discussed in the text as possible non-trivial fixed point. 

Fig. 3. Phase diagrams of U(l) systems with (a) naive (n/ = 2), (b) SLAC, (c) the 
Weyl, (d) mirror fermion (rif =2) actions. Dashed lines in (d) is extrapolation of the 
formulae ( |3.15|) to FI region. Points A in (a) and (6) are discussed in the text as possible 
non-trivial fixed points. 

Fig. 4. Phase diagrams of SU{2) systems with (a) naive (jif = 2), (6) SLAC, (c) the 
Weyl, (d) mirror fermion actions. Dashed lines is extrapolation of the formulae (|3.16|) to 
FI region. 

Fig. 5. Diagrams contributed to the fermion condensate ( |5.6| ) to the order (</>); (a) is 
the ladder diagrams. Solid circles stand for (</>), crosses for site n in ( |5.5[ ); other notations 
as in Fig. 1. 

Fig. 6. Fermion condensate [function C(y) in ( |5.6| )1 along the PM-FM critical line for 
Z2 system with SLAC fermions. 

Fig. 7. The same as in Fig. 6, but for U(l) systems with (a) naive, (6) SLAC fermions. 
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Fig. 2(c) 




Fig. 2(d) 
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